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ABSTRACT 


It  is  shown  that  a  one-way  3D  parallel/sequential  acceptor 
cannot  accept  the  class  of  3D  binary  arrays  in  which  the  set 
of  l's  is  connected,  unlike  the  situation  in  the  2D  case;  but 
it  can  do  so  if  multiple  passes  are  allowed. 
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0.  Introduction 


In  recent  years,  there  have  arisen  many  requirements  for 
three-dimensional  (3D)  data  processing  with  advances  in  com¬ 
puter  tomography  (CT).  Some  topological  properties  of  3D 
digital  pictures  have  been  discussed  in  a  series  of  papers 
by  Rosenfeld  and  Morgenthaler  [l]-[5],  as  well  as  others.  A 
3D  picture  can  be  represented  by  a  3D  array  of  volume  elements 
(voxels  for  short).  In  a  binary-valued  3D  picture  (each  voxel 
is  0  or  1),  it  is  easy  to  define  connectedness,  and  objects  and 
cavities  are  then  defined  as  the  equivalence  classes  of  the 
connectedness  relation.  These  correspond  to  two-dimensional 
(2D)  objects  and  holes,  respectively.  Moreover,  in  the  3D 
case,  there  also  exist  3D  holes  whose  properties  are  quite  dif¬ 
ferent  from  those  of  2D  holes. 

In  this  paper,  we  consider  the  problem  of  recognition  of 
the  above-mentioned  properties  of  3D  pictures.  First,  we  pro¬ 
pose  algorithms  which  for  every  (binary)  3D  digital  picture 
compute  the  numbers  of  objects,  cavities  and  holes.  These 
algorithms  are  performed  in  one  pass;  they  are  3D  versions 
of  the  algorithm  for  2D  digital  pictures  which  was  given 

by  Selkow  [6] .  Since  we  usually  receive  a  series  of  2D  pic- 

A 

tures  as  the  output  of  a  CT  scanner,  it  is  of  interest  that 
these  algorithms  scan  such  output  plane  by  plane,  from  the  top 
to  the  bottom,  and  stop  when  all  of  them  have  been  scanned. 
However,  the  algorithms  violate  the  condition  of  finiteness 
of  the  symbols  used.  From  the  point  of  view  of  the  theory  of 
languages  and  their  acceptance,  this  condition  of  finiteness 


is  a  fundamental  one.  In  the  second  part  of  this  paper ,  we 
discuss  this  problem  and  show  that  any  one-way  3D  parallel/ 
sequential  acceptor  cannot  accept  the  set  of  3D  pictures 
such  that  the  set  of  l's  is  connected.  This  is  interesting 
because  it  is  known  that  a  one-way  2D  parallel/sequential 
acceptor  can  do  this  acceptance  for  the  2D  case  (e.g.,  see 
Theorem  7.4.6  in  [7]).  In  the  third  part  of  this  paper, 
we  extend  the  concept  of  a  3D  parallel/sequential  acceptor 
to  allow  multiple  passes.  This  concept  of  multi-pass  is  the 
same  as  that  introduced  in  [8].  Then,  making  use  of  a  shrink¬ 
ing  technique  used  in  [9]  we  show  that  a  multi-pass  one-way 
3D  parallel/sequential  acceptor  can  accept  3D  connectedness. 

We  assume  that  the  readers  are  familiar  with  the  basic 
concepts  of  picture  languages. 


1.  Preliminaries 


Let  2  be  a  3D  array  of  lattice  points,  which  we  assume  to  be 
£xm*n,  i.e.,  2  =  {(i,j,k)|  l<i<£,  l<j5m,  l^k^n}.  A  3D  digital 
picture  f  is  a  mapping  from  2  to  {0,1},  i.e.,  f:2-*{0,l}.  Each 
point  (i,j,k)  is  called  a  voxel.  Here,  we  assume  that  any  3D 
digital  picture  is  surrounded  by  the  special  symbol  #  (called  the 
blank  symbol).  That  is,  let  B  be  the  set  {  (i , j ,k) |  (i=0  &  0Sj5m+l 
&  Osksn+l)  v  (i=£+l  &  Os j  £m+l  &  Osksn+l)  v  (Osi<£+l  &  j=0  &  Osksn+1) v 
(Osise+1  &  j=m+l  &  Osksn+1)  v  (0sis«+l  &  Osjsm+l  &  k=0)  v  (Osis^+l  & 
Osjsm+l  &  k=n+l)}.  Then  every  latice  point  of  B  is  occupied  by 
the  symbol  #.  The  set  B  is  called  the  border  of  2.  Usually, 
the  subset  of  l's  of  z,  i.e.,  the  set  {(i,j,k)|f(i,j,k)=l),  is 
called  S,  and  its  complement  is  called  S.  For  every  pair  of  points 
X=(x1,x2,x3)  and  Y= (y^ ,y2 ,y3 ) ,  X  and  Y  are  6-adjacent  if  |x1~y1|+ 

lx2-y2  l  +  lx3_y3  l=1;  y  and  Y  are  2 6-adjacent  if  max  (  |  x-^-y-J  ,|  x2-y2  |  , 

|x  -y  | )=1.  If  points  P  and  Q  are  6-adjacent  (26-adjacent) ,  then 
P  is  called  a  6-neighbor  (26-neighbor)  of  Q.  To  avoid  ambiguous 
situations  we  assume  that  opposite  types  of  adjacency  are  used 
for  S  and  S.  A  6-path  (26-path)  n  is  a  sequence  of  points,  tt  = 
p  ,Pl, . . . ,pm,  where,  for  all  i  such  that  lsism,  p^  is  a  6-neighbor 
(26-neighbor)  of  p^_^.  Any  two  points  P,Q  of  S  are  called 
connected  in  S  if  there  exists  a  path  P=Pg , . . . »Pm“Q  from  P  to  Q, 
where  p^S.  Evidently,  "connected"  is  an  equivalence  relation. 

This  relation  partitions  S  into  equivalence  classes.  These  clas¬ 
ses  are  called  the  connected  components  of  S.  In  the  same  way. 


we  may  define  connectedness  in  S  and  the  connected  components 
of  S.  A  connected  component  of  S  is  called  an  object  of  S. 
Clearly,  exactly  one  component  of  S  contains  the  #'s.  This  com¬ 
ponent  is  called  the  background  of  S;  all  other  components  of 
S  are  called  cavities  of  S. 

Even  in  ordinary  topology  it  is  difficult  to  characterize 
holes.  A  hole  may  be  thought  of  as  a  property  of  a  boundary 
surface  which  makes  it  topologically  equivalent  to  a  torus. 

In  another  approach,  an  object  is  defined  to  have  no  hole  if 
every  simple  closed  curve  in  the  object  is  continuously  deform¬ 
able  within  the  object  to  a  single  point.  We  see  from  these 
remarks  that  the  concept  of  a  hole  is  different  from  those  of 
objects  and  cavities;  we  cannot  point  to  or  label  the  points  which 
constitute  a  hole.  Indeed,  the  points  of  the  objects  and  cavi¬ 
ties  cover  the  space,  but  a  hole  is  a  property  of  those  collec¬ 
tions  of  points.  Thus,  when  considering  an  object  (and  its  cavi¬ 
ties)  we  shall  here  try  only  to  understand  what  is  meant  by  the 
number  of  holes  in  the  object,  and  not  what  is  meant  by  a  hole. 

The  genus  G (S )  of  a  set  S  in  a  3D  digital  picture  is  defined 
as  the  number  of  objects  in  S  (0(S))  plus  the  number  of  cavities 
in  S  (C(S) )  minus  the  number  of  holes  in  S  (H(S)).  As  already 
mentioned  the  definition  of  holes  is  not  simple,  and  in  parti¬ 
cular  holes  cannot  be  labelled  to  facilitate  counting  them. 

But  since  this  can  be  done  with  objects  and  cavities,  the  defi¬ 
nition  of  the  genus  defines  the  number  of  holes  in  S,  and  con¬ 


versely. 


In  [4] ,  Morgenthaler  has  given  a  method  of  computing 


G ( S)  directly  from  the  local  patterns  in  S: 

1)  When  26-adjacency  is  used  for  S, 

G26(s)“*i-*2+*3"*4+*5“*6+*7“V 

where 

4>1=#  (2) , 

4>2=#(3)+#(4)+#(5), 

4>3=#(6)+#(7)+#(8), 

<J>4=#(9)+#(10)+#(11)  +  #(12)+#(13)+#(14), 

4>5=#(15)  +  #(16)+#(17), 

4>6=#  (18)+#(19}+#(20), 

<J>7=#  (21 )  , 

<J>g  =  #(22), 

and  by  # (n)  we  mean  the  number  of  times  the  configuration 
n  of  Fig.  1  occurs  in  the  picture  S  (in  all  orientations) . 

2)  When  6-adjacency  is  used  for  S, 

Gs  (S)  =ip1-^2+4;3-ip4  r 
where 

^=#(2)  , 
i^2=#  (3)  , 
ip3=#  (9)  , 

ti>4=#(22)  , 

The  22  patterns  used  in  these  definitions  are  drawn  in  Figure  1. 
Morgenthaler  has  also  shown  that  G_c (S) -Gc (S) =1  and  Gc  (S) -Goc (S) =1 


.  The  number  of  objects,  cavities  and  holes 


In  this  section,  we  consider  algorithms  that  compute 
the  number  of  objects,  cavities,  and  holes  in  a  3D  digital 
picture. 

A  scanner  is  an  2.*m  array  of  finite-state  automata 
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Each  automaton  a.,  is  defined  by  a  7-tuple  <Q . . , 6 . . ,a . . ,a . . , 

13  13  13  13  13 

X..,&..,b..>  ,  where  Q. .,  the  set  of  states,  is  a  set  of  inte- 
x3  13  13  13 

gers,  and  the  next  state  function,  has  the  following  form: 

a .  . 

13 

6 .  .  s  n  Q  ...  *{0,1}  -»•  Q.  . 

13  k=1  aii(k)  13 

where  a^  is  an  integer,  and  cu  j : {1, . . . ,a^j}  -*•  1*1  is  a  one-to- 
one  function  which  enumerates  the  next  state  neighborhood  of 

0 i  . .  the  output  function,  has  the  following  form: 

b.  . 

13 

a.  . n  qr  m  x{o,i)  -  1, 

13  k=i  eij (k) 

where  b^j  is  an  integer,  and  j  :  { 1 , .  .  .  ,b^  ^  }  -*■  1*1  is  a  one-to- 
one  function  which  enumerates  the  output  neighborhood  of  ^ . 

We  will  use  ^  (t)  to  represent  the  state  of  0^  at  time  t. 
It  is  assumed  that  the  scanned  advances  by  one  plane  in  each 
unit  of  time  and  that  it  scans  the  first  plane  at  time  t=l. 
Thus,  the  input  to  scanner  element  at  time  t  is  f(i,j,t). 


Note  that  Q^j  is  not  a  finite  set  and  also  that  the  neigh¬ 
borhood  of  o. .  is  not  fixed. 

^  ] 

The  counter  C  monitors  the  output  of  each  element  of  the 

l  m 

scanner;  thus  C(t)=C(t-l)  +  E  Z  A..(t). 

i=l  j=l  13 

Now,  we  shall  describe  algorithms  for  computing  the 
numbers  of  objects,  cavities,  and  holes. 

1 :  Objects 

The  set  of  states  of  scanner  element  is  Q^j={x|x  is  an 
integer  and  |  x|<.(i+j-l)  (i+j-2) /2+j  }  .  Each  automaton  starts 
in  state  0  and  remains  in  that  state  as  long  as  0's  are  scanned 
When  a  voxel  containing  a  1  is  reached,  ck  j  will  assume  state 
(i+j-1) (i+j-2)/2+j.  As  tracks  a  string  of  l's,  an  exten¬ 
sion  of  the  component  of  S  is  sought,  i.e.,  two  automata  which 
are  actively  tracking  l's  are  tracking  the  same  object  if  they 
are  spatially  neighbors  or  if  they  are  in  the  same  state.  All 
automata  which  are  tracking  the  same  object  assume  the  same 
state  (the  state  of  the  automaton  having  the  smallest  state) . 

An  automaton  which  passes  a  lower  border  of  an  object  and  has 
been  in  a  state  k  enters  the  state  -k  for  one  period.  The  next 
time  it  goes  directly  to  state  0  unless  a  1  is  encountered.  If 
0^^  enters  the  state  -{ (i+j-1) (i+j-2) /2+j)  and  no  others  are  in 
the  state  (i+j-1) (i+j-2) /2+j ) ,  then  will  output  1,  i.e., 
one  object  has  been  scanned.  The  precise  definitions  of  6^ 
and  A. .  are  represented  as  follows: 


i)  When  6-adjacency  is  used  for  S, 
<5  i  j  (Q^  (t-1)  ,f  (t,i,  j) ) 


((i+j-1) (i+j-2)/2+j) -f (t,i, j)  if  Q.. (t-1)  £  0 


LINK ( Qi j (t-1) ) 
-LINK  (Qj^  j  (t-1)  ) 


if  (t-1)  >  0  a  f(t,i,j)  =1 
if  Qi j  (t-1)  >  0  a  f (t,i,j)  =0, 


1  if  Q^j (t)  =  -( (i+j-1)  (i+j-2)/2+j )  for  all  x,y  such  that 


0  otherwise, 


- (x+y-1) (x+y-2)/2+y  >  (i+j-1)  (i+j-2)/2+j , 

Qxy(t)*-Q..(t), 


where 


LINK (g ) 


=  min  2  ^  (xk+yR-l)  (xk+yR-2) /2+yR  |  (3xx)  .  . .  Ox^)  Oy^  . . .  Py^)  (Qx  (t-1) 

laksn  1  i 

g) A(Vj) ( 2s jsk) 

(  (Q  (t-1) >0) a ( (x . ,y . )  and  (x.  ,  ,y.  .)  are  6-adjacent) v (Q 

jYj  3  3  3~  3~  XjYj 

(t-1) =Q  v  (t-1)))))}. 

ii)  When  26-adjacency  is  used  for  S, 

5 .  • (Q._j (t-1) ,f  (t,i, j) ) 

1  J  J 


LINK( (i+j-1) (i+j-2) / 2+ j ) 
-LINK( (i+j-1) ( i+j-2) /2+ j ) 


is  the  same  as  in  the  case  (1) , 
where 


if  Q^j (t-1)  s  0  a  f (t , i  ,  j )  = 

if  f (t,i, j)  =  1 

if  Q^j  (t-1)  >  0  a  f  ( t ,  i  ,  j  )  =0 


LINK (g) 


2  {  (xk+yk_1)  (V"yk~2)/2+yk  1  (3xl)  •••  (3xk_i>  (3v1)  •••  <3yk_1)  ((x1+y1-l) 

•(xl+yl~2)  /2+yi 

=g) a ( V j ) (2sjsk)  ( (Q  (t-l)>0) a ( (x . ,y . )  and  (x.  .,y.  .)  are 

xjYj  j  j  D D 

26-adjacent) v (Q  (t-l)=Q  (t-1) ))))}. 

jYj  Xj-lYj-l 

2 :  Cavities 

Since  all  components  of  S  except  the  background  component  ar 
cavities  of  S,  the  algorithm  for  counting  objects  of  S  can  be 
used  for  counting  cavities  of  S  by  interchanging  the  roles  of  1 
and  0.  In  this  case,  the  initial  value  of  the  counter  C  must  be 
to  remove  the  background  component  from  the  set  of  cavities  of  S. 

3 :  Genus 

For  any  S,  every  2*2*2  local  pattern  in  Fig.  1  is  easily 
counted  by  our  computational  model.  Thus  the  algorithm  computing 
the  genus  of  S  is  easily  defined. 

4 :  Holes 

Finally  from  algorithms  1-3,  we  can  construct  the  algorithm 
for  counting  the  number  of  holes  in  S,  since  H (S) =0 (S) +C (S) -G (S) . 


deterministic  one-way  PSA 

In  Section  2,  we  have  given  algorithms  which  for  every  3D 
digital  picture  count  the  number  of  objects,  cavities,  and  holes. 
However,  as  noted  in  Section  2  the  used  in  the  algorithms 
needs  a  set  of  states  which  is  not  finite.  Further,  its  neighbor¬ 
hood  is  not  fixed.  In  language  theory,  finiteness  is  required  in 
the  automaton  used.  From  the  point  of  view  of  language  theory, 
we  prove  here  the  unacceptability  of  the  connectedness  of  3D 
pictures  by  deterministic  one-way  parallel/sequential  acceptors. 
Parallel/sequential  acceptors  of  2D  pictures  are  well-known  (e.g., 
see  [7]).  It  is  also  known  that  a  deterministic  one-way  non¬ 
writing  parallel/sequential  acceptor  can  determine  whether  or  not 
the  set  of  x's  in  a  rectangular  array  Z  of  x’s  and  y's  is  connected. 
In  contrast  with  this  result,  our  new  theorem  for  3D  seems  to  bf 
interesting . 

A  parallel/sequential  acceptor  of  3D  pictures  (for  brevity: 

3DPSA)  is  analogous  to  that  of  2D  pictures.  It  is  defined  as 
follows:  A  parallel/sequential  acceptor  is  a  9-tuple  A=(Q,qQ,QA, 

#, V, #  , #b, 6 ,y ) ,  where 

Q  is  a  finite  nonempty  set  of  states 

qQ  (  Q  is  the  initial  state 

Qa  c  q  is  the  set  of  accepting  states, 

#  €  Q  is  the  blank  symbol, 

V  is  a  finite  nonempty  set  of  symbols  called  the  tape  vocabulary, 
(L  and  #,  are  blank  symbols  in  V, 


6:  Q"xv  2“  ’  is  the  state  transition  function  and 
p:  Qxv  -*■  2 ^  is  the  move  function. 

The  operation  of  A  on  a  3D  array  E  can  be  described  as  follows: 
A  consists  of  a  2D  array  of  cells 


C11  ••• 

Clm 

C21  •** 

c2m 

C  p  1  •••  C 

^1  ^m 

whose  lengths  are  the  numbers  of  columns  and  rows  of  E,  respective¬ 
ly,  together  with  special  "cells" 


l+l  0 


that  are  regarded  as  permanently  in  the  #  symbol.  E  has  a  row 
of  #  ’ s  just  above  its  top  array  and  a  row  of  #b’s  just  below 
its  bottom  array.  Initially,  A  is  on  the  top  array  of  E  with  every 


c00  C01  *  *  * 

C0 

m+1 

c10 

C1 

m+1 

C20 

C2 

m+1 

• 

• 

• 

• 

'S-m+1 


'£+1  m+1 


cell  in  the  state  qQ.  At  any  given  step,  each  cell  c^j  reads  the 
symbol  v  in  its  position,  senses  the  states  *  *  *  '^9  ci-l'j- 

ci-l,j'  ci-l,j+l'  ci, j-1 '  Ci,  j '  Ci, j+1'  Ci+1 , j-1 '  ci+l,j+l' 

and  can  go  into  any  new  state  q’  and  write  any  new  symbol  v'  such 

that  (q' ,v' ) €6 (q1,q2f • • • fqg) •  The  move  function  depends  only  on 


the  (new)  state  and  symbol  read  by  the  distinguished  cell  c^. 

0€p  means  that  A  can  stay  where  it  is;  l€p  means  that  A  can  move 
down,  -l€p  means  it  can  move  up. 

It  is  required  that  M(q,#t)=l  and  U (q, #b) =  -1  for  all  q.  It 
is  understood  that  #t  and  can  never  be  rewritten.  If  c^  ever 
enters  a  state  in  QA,  we  say  that  A  has  accepted  E. 

If  a  PSA  does  not  move  up,  A  is  called  one-way. 

Now,  we  prove  the  following  theorems.  Let  C  be  the  set  of 
3D  arrays  consisting  of  O's  and  l's  in  which  the  l's  are  connected. 

Theorem  2 . 1  C  is  accepted  by  a  nondeterministic  one-way  3DPSA. 
Proof:  Let  us  consider  isometric  3D  array  grammars,  which  are 

three-dimenisonal  analogs  of  those  in  the  2D  case  (see  [7]). 

Notice  that  C  is  generated  by  an  isometric  3D  monotonic  array  gram¬ 
mar  G.  Here,  we  can  assume  that  the  starting  symbol  p  of  G  appears 
only  once  during  the  applications  of  rewriting  rules. 

Now,  let  A  be  a  nondeterministic  one-way  3DPSA  working  as 
follows : 

1)  When  reads  a  symbol,  the  state  changes  into  one  corres¬ 
ponding  to  a  non- terminal  symbol  of  G.  This  change  is  done 
nondeterministically . 

2)  Illegal  guessing  in  1)  causes  o^j  to  go  into  a  dead  state. 

3)  When  a^j  goes  into  the  state  corresponding  to  the  starting 
symbol  p,  this  state  moves  to  o^*  In  this  case,  A  never 
goes  downward  until  this  move  is  completed. 

4)  For  two  states  p,  a-.  goes  into  the  dead  state. 


5)  When  A  reaches  the  bottom  plane  and  finishes  its  transitions, 
°11  9oes  into  the  accepting  state  if  it  has  memorized  the 
only  p. 

From  the  above  construction,  we  have  the  theorem.  || 

Theorem  2.2  C  is  not  accepted  by  any  deterministic  one-way 
3DPSA. 

Proof:  Let  us  consider  a  square  array  of  side  length  4Z  as  shown 
in  Fig.  2.  The  left  and  right  halves  are  called  the  L-part  and 
R-part,  respectively.  Cells  of  both  parts  are  occupied  by  0's 
and  l's  as  shown  in  Fig.  3.  That  is,  every  even  row  and  every 
even  column  are  occupied  by  l's,  and  all  other  cells  by  0.  This 
array  is  put  in  some  plane  of  Z.  In  the  arrays  other  than  this 
plane,  the  1-cells  are  connected  as  shown  in  Fig.  4.  Here,  n^ 
is  the  name  of  a  cell  in  the  L-part  and  m^  is  the  name  of  a  cell  in 
the  R-part.  Thus,  the  cells  ni  and  n^+1  in  the  same  parentheses 
are  connected  in  the  cells  between  the  top  array  and  pth  array, 
and  the  same  for  the  m^  and  nu+1  cells. 

Further,  the  line  indicated  in  Fig.  4  shows  that  the  n^  and 
itk  cells  are  connected  in  the  cells  under  the  pth  array.  This  con¬ 
nection  is  always  possible  by  considering  a  very  large  p  and  suff- 
ciently  high  Z . 

Now,  the  number  of  configurations  of  states  of  A  at  the  time 

4  8,  x  4  8, 

it  leaves  the  pth  array  is  | Q |  .  Also,  the  number  of  connec¬ 

tions  of  cells  in  the  L-part  (that  is,  the  number  of  parenthesized 


9^  2  r  21  1 

1-cells)  is  2JL  2  2)1  -2  2  ...2  2  =  - 51-  .  But  we  have 

2^ 

2  2  z 

| q | 16 £  ^  -  for  sufficiently  large  i.  Hence,  there  are  at 

2* 

least  two  different  parenthesizings  which  yield  the  same  con¬ 
figuration  of  states  of  A  at  the  time  it  leaves  the  pth  array. 

Let  these  two  be  as  follows: 

...  (n^ ,  n0 )  (n^  »n^ )  ... 

...  (m^  f  m^  ^  f  m4  ^  ... 

and 

...  (n^ ,n^ )  (  )  ... 

...  (m^ f m^ )  (m3 > )  ... 

Let  Fig.  5  and  Fig.  6  represent  these  two  3D  arrays,  where  the 
connection  relations  are  in  the  cells  between  the  (p+l)st  and  bot 
tom  arrays.  The  lines  in  Fig.  5  show  that  all  1-cells  are  con¬ 
nected.  This  is  possible  by  considering  a  suitable  connection 
level. 

Suppose  that  a  deterministic  one-way  3DPSA  A  accepts  the 
connected  3D  arrays.  Then  A  accepts  the  connected  3D  array  shown 
in  Fig.  5.  But  then,  A  must  also  accept  the  non-connected  3D 
array  shown  in  Fig.  6,  since  it  gives  rise  to  the  same  configura¬ 
tion  of  states  of  A  at  the  time  it  leaves  the  pth  array.  This 
is  a  contradiction.  Therefore;  we  have  the  theorem.  // 


4.  Acceptability  of  the  connectedness  of  3D  pictures  by  a  multi- 


pass  one-way  PSA 

In  the  previous  section,  we  showed  that  any  deterministic 
one-way  PSA  cannot  accept  the  connectedness  of  3D  pictures.  Here, 
we  modify  the  one-way  PSA  so  it  can  make  repeated  passes  over  a 
given  3D  picture.  That  is,  the  new  acceptor  A  acts  as  follows: 

1)  First,  A  works  as  the  usual  one-way  PSA.  This  is  the  first 
pass . 

2)  When  A  reaches  the  array  below  the  bottom,  A  begins  to  work 
again  from  the  top  array.  This  is  the  second  pass. 

3)  A  repeats  the  behavior  2). 

4)  If  c^  ever  enters  a  state  in  QA  during  the  repeated  behavior 
of  3),  we  say  that  A  accepts  E. 

For  this  multi-pass  3DPSA,  we  have  the  following  theorem: 
Theorem  3.1  A  multi-pass  deterministic  one-way  3 DPS A  can  accept 
the  connectedness  of  3D  pictures. 

Proof:  Let  u(t)=l  if  t>0,  u{t)=0  if  t<0.  We  consider  3D  pictures 

consisting  of  0's  and  l's.  In  Fig.  7  a,b,c ,d,e, f ,g,h  are  the 
voxels  belonging  to  a  2x2x2  window. 

Then  F(b),  the  transformed  value  of  b,  will  depend  on  the 
values  of  the  elements  belonging  to  the  three  planes  that  meet  at 
B. 

For  6-connectivity,  we  define 
F (b) =u (u (a+b-1 ) +u (b+c-1 ) +u (a+d+c-2 ) ) +u (f +b-l ) +u (c+g+f-2 )+ 
u (a+e+f-2 ) ) . 

Similarly,  for  26-connectivity,  we  define 
F (b) =u (u (a+b+g-1 ) +u (b+c+e-1 }+u(b+f+d-l)+u (a+c+f-1 ) + 


Then  the  following  theorem  is  well-known  [9] :  When  F  is  applied 
repeatedly  in  parallel,  a  single  object  shrinks  down  to  a  single 
1-voxel  in  a  finite  number  of  steps. 

By  making  use  of  this  result,  we  can  prove  the  theorem. 

This  is  done  as  follows: 

In  the  first  pass  of  a  multi-pass  deterministic  one-way  PSA  A, 
the  values  of  a,c,d,e,f,g,h,  are  written  at  b  as  the  output  of 
A.  This  is  done  for  all  voxels. 

In  the  second  pass,  A  writes  the  result  obtained  by  the 
first  application  of  the  function  F  to  E.  Generally,  in  the 
(n+l)st  pass  A  writes  the  result  obtained  by  the  nth  application 
of  F  to  E. 

When  A  recognizes  a  single  1-voxel,  c^1  enters  the  accepting 
state.  Therefore,  we  have  the  theorem.  // 

By  a  similar  technique,  A  can  accept  E  such  that  it  has  a 
cavity.  Furthermore,  by  making  use  of  a  counter  (bounded  to  the 
number  of  voxels  in  E )  a  multi-pass  deterministic  one-way  PSA 
A  can  count  the  number  of  objects  and  cavities  in  E.  Thus,  we 
know  from  the  equation  given  earlier  that  a  multi-pass  determini¬ 
stic  one-way  PSA  can  count  the  number  of  holes  in  E. 
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